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1 Introduction 



String theory has long been known to be a generahzation of gauge theory due to the presence 
of a massless spin one state in the open string spectrum. Since all of the massive states of 
the theory have masses proportional to \/2t^Tq — 1/ \fa', the string theory goes over to the 
gauge theory in the infinite tension limit a' — )■ 0. This, together with the fact that closed 
strings have a massless spin two state, has inspired the discovery of many deep connections 
between string theory on the one hand and gauge theory coupled to gravity on the other. 
The AdS/CFT correspondence [1] which asserts the equivalence of = 4 supersymmetric 
gauge theory to type IIB superstring theory on an AdSsxS^ space time manifold is one 
of the most spectacular of these. Although motivated by the physical properties of open 
string theory with a' > 0, the final conclusion is reached by taking a' — )■ 0. The finiteness 
(conformal invariance) of the J\f = 4 theory plays a key role in justifying the a' ^ limit. 

The corresponding hypothesis for an asymptotically free gauge theory like the gluonic 
sector of QCD is more obscure. However, there is little doubt that string theory can offer 
important insights into some aspects of QCD. In this article we launch a critical analysis 
of the possibility that 't Hooft's — )■ oo limit [2] of QCD might be usefully analyzed 
by replacing it with the sum of open string planar diagrams, keeping a' > 0. There are 
several reasons to hope this helps. First, the organization of multiloop string diagrams is 
dramatically simpler than the corresponding gauge theory diagrams: there is only one planar 
open string diagram at each loop order, whereas the number of planar gauge theory diagrams 
grows exponentially with order. Secondly, the loop open string planar diagrams describe 
the evolution of a "bare" worldsheet which becomes "dressed" with the inclusion of planar 
loops. This provides a very natural setting for the description of a confining fiux tube which 
may survive the hmit a' — > 0. Finally there is the long held expectation that the ultraviolet 
behavior of gauge theory diagrams will be mitigated by the "stringiness" associated with 
finite a', making the latter better defined. 

We focus on this last point in this article. It is not so much the ultraviolet divergences 
themselves that concern us here-after all those can be absorbed in coupling renormalization 
in gauge theories. Rather, it is the extreme care that must be taken in gauge theory to 
preserve gauge and Lorcntz invariance in the finite part that remains after renormalization. 
Order by order in perturbation theory, this can be accomplished by employing a suitable 
regularization, the most popular of which is dimensional regularization. But as soon as one 
aims to extract the consequences of summing all the planar diagrams, especially if one must 
rely on numerical methods, the soundness of dimensional regularization becomes somewhat 
questionable. For example the powerful conclusions derived from lattice gauge theory would 
be much less convincing if they relied on an unphysical regularization such as "analytic 
continuation" of the dimension of spacctimc. It is desirable to have a digitization scheme 
which can be relied on to give the correct physical results without such an artifice. Thirty 
five years ago Giles and one of us [3] (GT) proposed a digitization of the sum of planar 
open string diagrams based on lightcone quantization [4, 5] in its path history formulation 
[6] . In the present article we set out to assess the reliability of this specific lattice model for 
perturbative calculations. 
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Let us briefly review the GT proposal, in order to set the stage for the rest of the 
paper. In lightcone quantization of the bosonic string [5], one takes x+ = (x° + x^)/^/2 as 
the quantum evolution parameter and labels points on a string by a parameter o defined 
so that P"*" = (P° + P^)l\/2 is uniformly distributed on the string. In effect these two 
choices eliminate and x" as dynamical variables, leaving only the D — 2 transverse 
coordinates x{a) as quantum operators. Mandelstam worked out the path history form of 
this quantization [6] using imaginary time r = Then the propagator for a free string is 
simply the path integral over the x{a^ t) where < a < P+ and < r < T parameterize 
a rectangular region of dimensions P+ x T. The path integrand is simply , with S the 
lightcone Euclidean action 



In this language a general open string planar diagram is calculated by integrating this same 
integrand over a worldsheet with several slits of variable length and location as depicted in 



T 

Figure 1: Typical open string planar diagram on the lightcone worldsheet. 
This one is a seven loop 5 string function 

The GT proposal is simply to digitize Mandclstam's interacting string diagrams by defin- 
ing a rectangular M x N grid with T = (A^ + l)a and P"*" = MoTq. Then the integration 
variables a;((T, r) a?/ and the lattice action is simply 




(1) 



Figl. 



P+ 
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A quick look at the lattice corresponding to a multiloop open string diagram Fig. 2, shows 
that a slit is nothing but a row of missing spatial bonds (links), and summing over all 
planar diagrams is simply summing over all patterns of missing spatial bonds. One can 



Figure 2: Multiloop lattice worldsheet: each loop is a row of missing links. 

easily incorporate the sum over missing bond patterns into the lattice sum over histories by 
introducing an Ising-like variable = 0, 1 on each spatial link. Then the worldsheet action 
describing the sum of all open string planar diagrams is simply 

ij 

The purpose of the In g term is to insert a factor of g at the beginning and end of each row 
of missing bonds. Then, in addition to integrating each from — oo to +oo, one sums each 

over the values and 1 to obtain the sum over all planar diagrams. 

If no further adjustments were necessary, we could study this system numerically, e.g. 
through Monte Carlo simulation. By analyzing the large behavior of the path integral, one 
could read off spectral information by identifying exponential behaviors e~"^^(^)^, where 
Ex{M) are the eigenvalues of P~ . On general grounds we should expect the large M behavior 

Ex{M) - aM + /3+^ + --- (4) 

The aM term is a bulk worldsheet effect and the /3 term is associated with boundaries. The 
coefficients a, /3 depend on the details of the lattice model and violate the requirement that 
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2p-p+_p2 _ -g Lorentz invariant. Fortunately there are two geometrical counterterms 

naturally associated with the worldsheet path integral: one proportional to the area of the 
worldsheet and another proportional to the length of worldsheet boundaries. Thus the 
Lorentz violating terms noted above can always be cancelled. As noted in [3] the area term 
is dynamically inconsequential for the sum of diagrams because the slits representing loops 
have zero area, so that the area term is identical for all diagrams contributing to the same 
process. The boundary term depends on the number and lengths of the slits. It can be 
thought of as an energy cost assigned to the disappearance of a bond. In the Ising spin 
description of the sum over loops it is represented by a term {D — 2)B — S/) added 

to S. The boundary term is absent for closed strings, but is required for a Lorentz invariant 
free open string spectrum. This is because the Gaussian lattice worldsheet path integral 
implies an open string zero point energy 



as M — )■ oo, where G is Catalan's constant. Cancelling the second term requires B ^ Bq = 
(1/2) sinh~^(l) for g = Q. Since the l3 term in the energy will in general depend on g, we can't 
know the value of B at finite g a priori, so in practice it must be left as a free counterterm 
parameter to be determined by requiring that the final answer be consistent with Lorentz 
invariance. 

More generally, we can regard B as a free parameter of the lattice model, which we don't 
necessarily have to insist is Lorentz covariant. Taking B large enough lifts the P~ of any 
multi-open string intermediate state above all the g = closed string energy eigenstates that 
survive the continuum limit^. Then the severe infrared divergences in loop diagrams caused 
by the open string tachyons are removed, making the multi-loop expansion well defined, 
albeit with a loss of Lorentz covariance. For numerical studies we should therefore calculate 
for general B large enough to remove instabilities; and only at the end of the calculation 
would we scan for a value of B which restores Lorentz covariance, if possible. 

But there is no a priori guarantee that these two counterterms can remove all Lorentz 
violating artifacts from this lattice worldsheet construction. For instance, the authors of 
[7] developed a lightcone worldsheet formalism which mapped the planar diagrams of gauge 
theories in lightcone gauge to a worldsheet system with exactly the features of the open string 
planar diagrams just described. String coordinates were employed, but their dynamics were 
topological in the sense that all but a single zero mode decoupled in the path history sum. 
In this way the worldsheet diagrams were designed to yield precisely the "bare" Feynman 
diagrams of the gauge theory. Digitization of this worldsheet in the manner of GT amounted 
to a scheme for cutting off the UV and IR divergences of these Feynman diagrams. To test 

''in effect with B this large one imposes confinenient on the free theory! The crucial issue is then whether 
or not confinement survives as B is reduced to a value that restores Lorentz invariance. 



M-l 



aP-{M) 



2) ^ sinh-^ sin WT7 + {D - 2)B 



m=l 
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the reliability of this cutoff, the authors of [8] calculated the one loop diagrams contributing 
to the scattering of glue by glue regulated by the GT worldsheet lattice. They found that the 
artificial divergences associated with lightcone quantization could indeed be absorbed in the 
bulk or boundary worldsheet counterterms. Further the divergences associated with charge 
renormalization (asymptotic freedom) had the correct coefficients. Unfortunately, this was 
not the end of the story. There remained gauge violating terms that could only be cancelled 
by: (1) a divergent gluon self mass, (2) a finite wave function renormalization, (3) a finite 
adjustment to the three gluon function, and (4) a finite constant adjustment to the 4 gluon 
function. Indeed these are precisely the adjustments generally required when a physical 
cutoff is employed in loop calculations [9]. All of these adjustments, being polynomials in 
the external momenta, are consistent with locality. The bottom line is that the lightcone 
worldsheet lattice as a regulator of gauge theory diagrams is no better than other physical 
cutoffs. Unfortunately this means that we should expect the necessary counterterms to 
proliferate with the inclusion of multiloop diagrams. 

In this paper we begin to explore whether the GT lattice does a reliable (or at least better) 
job regulating open string planar diagrams than it does with field theory planar diagrams. 
Does keeping a' > control the proliferation of counterterms? Is it possible that the bulk 
and boundary counterterms will be sufficient by themselves? We start with the simplest 
self energy diagram: the one loop correction to the closed string propagators. This process 
involves only a single intermediate open string state and hence has the singularity structure 
of a tree diagram. Since the corrections to the open string propagator involve complications 
associated with the multi-string intermediate sates and worldsheet boTindarics. wc choose to 
defer the open string analysis to a subsequent paper, and we restrict our attention here to 
the corrections to the closed string propagator. 

Because of the tree structure of the closed string self-energy diagrams, the sum over 
K, the number of time steps that the intermediate open string exists, converges for both 
the closed string tachyon and graviton, even with B = Bq. Nonetheless the calculations of 
Sections 3 and 5 (with B = Bq) will establish Lorentz violations in both the tachyon and 
graviton. Fortunately, these Lorentz violations disappear if the calculations are done holding 
B > Bq, taking B ^ Bq only at the end of the calculation. In any case, this is the only 
way to make sense of multi-loop diagrams, so we don't think its necessity at one loop is a 
drawback. 

The article is organized as follows. In Section 2 we obtain explicit formulas for the one 
loop energy shifts of the closed string ground state (tachyon) and the closed string graviton 
state. Each term in the K summand involves determinants of M x M overlap matrices. 
These formulas are evaluated and analyzed with the help of Mathematica in Section 3^. 
In Section 4 we extend the formulas to include open strings ending on Dp-branes, and 
similarly perform their numerical evaluation in Section 5. In this case, instead of an energy 
shift, the one loop correction gives the amplitude for a closed string scattering off the Dp- 
brane. Additional discussion is given in the concluding Section 6. Several appendices collect 

^For the interested reader, we provide the evaluation code and a sample of our analysis in a Mathematica 
file accompanying the source format of this article on the arXiv. 
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a number of technical results, including normal mode expansions, determinant formulas, 
and overlap matrices, that are useful for the discussion in the main text. We also include 
additional evidence for the robustness of our numerical results. 



2 Closed String Self-Energy 

Before calculating the self-energy diagram on the lattice, we briefly recall the known expres- 
sion for the continuum self-energy. In the language of conformal field theory we need the 
amplitude with two closed string vertex operators, say at 1 and oo, on the complex plane 
from which a disk of radius g < 1 has been excised. Consulting for example [10], we find the 
result for the closed string tachyon self-energy: 



which is obviously seriously divergent at q — 0: there are quadratic and logarithmic diver- 
gences in the q integration. To get further insight into the fate of these divergences, we do 
the path integral in lightcone parameterization, using the conformal transformation methods 
of [11, 12]. The (still seriously divergent) result in D spacetime dimensions is 



-AP- = C"P+ / dT 



j 

Jo 



27r 



P+sinh(7rrTo/P+) 



(D-2)/8 

(7) 



where T is the length of the slit on the hghtcone Mandelstam diagram, which is the total 
ix'^ over which the intermediate open string propagates. The factor of P'^ is just the result 
of integrating the a independent integrand over vertical location of the slit < o" < P^. 
When D — 26 this expression reduces to (6) which can be seen with the change of integration 
variable 



1 _ g-7rrTo/P+ 

1 + e-^^^o/-P+ ' 



(8) 



We can get a rough idea of what we should get from the lattice calculation by simply 
discretizing T — aK and P+ — qMTq. Then the discretization of (7) reads: 



-aAP- = (C"ro-2(aro)(26-^)/s)M5^ 



K=l 



2n 



Msmh{7rK/M) 



(D-2)/8 



(9) 



Of course discretizing the result of a continuum calculation is not the same as doing the 
discretized calculation from the beginning, but it at least can serve to guide the eye. For 
example one feature we immediately see from (9) is that the quadratic divergence seen in 
(6) is expected to be reflected in the lattice calculation as a linear term in M in P~, which 
can be absorbed in the bulk counterterm described in the introduction. We turn next to the 
actual lattice calculation of the self-energy. 
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M 



J K L 

Figure 3: Lattice worldsheet for the closed string self-energy. The dotted lines 
are identified. There are K — 1 missing links. 



We start with the expression for the summand of the one loop correction to the closed 
string propagator, depicted in Fig. 3. It is a product of factors, one for each of the D — 2 
transverse coordinates xf. In the following we will display only 1 of the factors, calling its 
coordinate x^: 



{N+l,{xf}\0,{x'}y"'"'''{K, J) 



closed/' 



dxfdxfiL, {xf}\0, {x'^jr^-^^'iK, {x^llO, {x^}r^- 



where J + K + L = N + 1, and Bq — (1/2) sinh~^ 1 is the counter term which removes the 
boundary contribution to the free open string P~ . Regarding the intermediate open string as 
a closed string with a row of missing links, we see that — 1 is the number of missing links, 
which we have taken to be the ones hnking site 1 to site M. The open and closed string free 
propagators are defined in Appendix C. The factors V, which are the corresponding propa- 
gators for vanishing values for the initial and final coordinates, are related to determinants 
of the worldsheet discretizcd Laplacian, and arc defined in Appendix B. 



Energy eigenvalues are determined by identifying exponential time (r 



IX 



in the closed string propagator e ^ ^ e 
and 



-a{N+i)p perturbation theory P' 



g-a(iV+l)p- ^ g-a(iV+l)p- [1 _ + ;l) AP- + • • •] 



) dependence 
= P,-+AP- 



(11) 



so we identify —aAP as the coefficient of + 1 in the one loop correction to the closed 
string propagator. The factor A" -|- 1 is associated with one of the sums over the creation 
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and destruction times of the intermediate open string. In practice the factor is removed by 
using time translation invariance to fix, say, the creation time, and summing only over the 

destruction time. 

If we desire to use this formalism to calculate the corrections to excited energy levels, it 
is necessary to deal with the fact that then the summand over K grows exponentially with 
K. A convenient approach is to add a term —Ke to the exponent, which can be thought of 
as adding e to i?o- Since all of the physically meaningful excited states have an excitation 
energy of order 1/M, as long as e > is independent of M this is enough to tame the 
exponential divergence. Actually, since Bq is a necessary counterterm which is expected to 
receive corrections in perturbation theory, we cannot know a priori what its value should be. 
All we know is B — (1/2) sinh~^ 1 + 0{g'^). Thus in the context of nonperturbative studies 
of the worldsheet path integral, it should be taken as a parameter to be tuned at the end of 
the calculation to ensure Lorentz covariance. 

For some states, such as the ground state and graviton state, the absence of lower energy 
open string states that couple means that the K sum converges without the need for e. 
However, as we shall see later, because the convergence is then only hke e~^/^, ultraviolet 
divergences become entangled with infrared divergences, which then leads to violations of 
Lorentz invariance since M oc P+. Introducing an e > prevents these artifacts from 
entering. Thus, we should keep e > even in those fortuitous cases where it isn't strictly 
necessary for convergence. 

2.1 Correction to the Closed String Ground Energy 

Since there is only one ground state, it is uniquely singled out by taking J, L ^ oo, in which 
case it suffices to simplify matters by setting = — 0. Then iW is simphfied to 

m m 
+ E(^- + COth KX:^-2J2 ^omQom^^] (12) 

In this equation AJ^ = 2 sinh^"*^ sin(m7r/2M) and A^^ = 2 sinh~^ sin(m7r/M) are the discrete 
time versions of the normal mode frequencies for the open and closed strings,respectively. 
Correspondingly, the g's are the normal mode coordinates of the open or closed string, defined 
in Appendix A. Note that when it makes no essential difference, we shall restrict M to be 
odd to keep the description of the closed string modes as simple as possible. 

We change integration variables to the closed string normal mode coordinates Qq, qcm, Qsm 
for both K and L. The Jacobian for this variable change is unity, and we can express the 
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^om in terms of the closed string modes as follows: 

?cm/2 

Q.om 



2^ 

M 



^{M-l)/2 rj 



m even 



m odd 



(13) 



where the overlap matrix Umm' is defined in Appendix D. We also need 

■(M-l)/2 



xm — xi — —2 




M 



E 

m'=l 



Qsm' sm 



w 



(14) 



Then 



Tn 



- f [(a^i. - )^ + (x^ - xf )^] 



To 
" 2 



J X 



+ Yl 'icm sinh A^(coth JA^ + coth iTA^) + ^ g^;^ sinh A^(coth LA^ + coth KX'J 



'sinh KXt, 



+ XI ^^^^ ^^^^ ■T'^m + X! 9sm sinh A^ coth LA^ 



+ E (5o^ + 9ot')sinhA;;COthi^A^-2 ^ ^on^^ot- 



sinh A° 



m odd 



m odd 

m'vr m"7r 



sinh XA° 



+ M ^ (5sm'5.m" + Qsm'Qsm") sin ^ sin 



M 



(15) 



Because of the equality qo2m = Qcm, the integration over qq and the Qcm precisely implements 
closure on these modes. This means that the result of those integrations is just the contri- 
bution of those modes to p^iosed^^ _|_ Consulting Appendix B for the various Z>'s, the 
outcome of the integration over go and the qcm can be written 



{N + l,{xf}\0,{x'}y^°'^'^ 



L AW'-{K-l)Bo 



(16) 



2 



Y,{(lfm) Sinh A^ coth JA^ + ^(^.m) sinh A^ coth LA 



+ E +?L'?W')(l^^ E t^mm'C/mm" sinh A^ coth XA^ 



VM2 



m odd 



m'TT 



m"Ti 



sm 



M 



m odd 



sinhA^ 
sinh /TA^ 



(17) 



and we remind the reader that there are D — 2 such factors. To isolate the shift in a specific 

energy level, we need to identify the exponential behavior in J, L as they approach infinity. 

For the ground state energy shift, it is sufficient to directly take J, L ^ oo with K fixed. 
Then the coth J, coth L ^ 1 and 



iCloscd 
sin 



(iV+1) ^ e 



-) 

27ry 



'ncloscd/ TX'ncloscd / r \ 



{M-l)/2 ,c 



27r/ 



(M-l)/4 (M-l)/2 

V2sinh 

m=l 

(M-l)/4 (^-l)/2 

JJ V2sinhA^ 

m=l 



Change integration variables to = Qsm^/ (^o/vr) sinh X^. And we have 
{N + l,{x^}\0, {x^yy^"'''^ 



m 



(18) 
(19) 



^(^o/Tr) si^IhA^ 



closed 



(7V+1) 



n 



llm=l,odd V -L e - ^^-^ 

where iH^" is defined below. It can be shown that 

M-l 



'sinh Aj^ 
sinh A^ 



cig^cig4e^^"-(^-^)^° (20) 



HsinhA^ = 2-(^-^)Wm 



m=l 



M-l 



sinh 2M sinh"^ 1 
sinh 2 sinh~^ 1 



2-(^-i)Wm 



sinh M sinh ^ 1 cosh M sinh ^ 1 



V2 



JJ^ sinhA^ 



2-(^-i)MsinhMsinh"^l 



(21) 
(22) 



m=l 



SO, restoring all D — 2 factors in the amplitude and summing only over K, the time interval 
spanned by the open string propagator, we infer 



-aAP, 



G.closed 



K 



(iV + l,(0)}|0,{0}) 

X)closed(7V _^ 1) 



closed' 



D-2 



- MY. 



eKj:Zzti^'k-^?.)/2-{K^i)Bo /cothMsinh-1 1 



-2KX° 



MV2 



nM-i r\~z — 
m=l,odd V J- e 

2 L^-^^- ^ ^ ^ v^iEhA^;;^EiA^ 

( 4 . 1 ,1 T^xo 4 . wl-K . m"v:\ 

1^ [/w C/mm" smh A„ coth irA^ + —sm-^sm^j 



D-2 



(23) 



m odd 
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M2 



Q^m'Qsm" TT TT sinhA^ 

/ ^ '-^ mm' mm" ' 



y - 

St" V^inhA^, sinhA^,, 



m odd 



sinh KX<L - 



= — TT 



TO' ?to" + Q.m'Q.m")^rn'm" 1m''lm"^rn'm" 

m',m" to' 

With the definitions on the last fine 

= det-V2(^^ ^) =det-V^(A + S)det-V2(A_S). 

The last equahty follows because the eigenvectors of 



/ 



(24) 



(25) 



A B 
B A 



can be taken to be of the 



form 



^ ^ ^ where v± is an eigenvector A± B. Finally we summarize the result 



-aAP, 



GjClosed 



MY. 

K=l 



Afn'm'' 



/ cothMsinh-^l \ e^i:^!^i(^^.-^^)/2-(^-i)go 
<^TO'm" _|_ 2 sin (mV/M) sin (m'V/M) 



det-^/2 



A B 
B A 



D-2 

(26) 



2 M sinh A^, sinh A 



Bm'm" 



TO odd 
Umm'Umm" sinh A^ 

M2 x/sinh A^, sinh A^^„ sinhTTA^ 



2 ^ ^toto" sinh A^ coth i^A^ 



-ysinh A^, sinh A^„ 



E 



(28) 



(^ ± -B)to'to" 



Srn/7 



2 sin (m'TT /M) sin (m'V/M) 



2 M ^ sinh A^, sinh A 



2 U U mm" sinh A^ 

^TO^dV^^TO'SinhA^,, 



tanh 



i^Al 



±1 



(29) 



2.2 Correction to the Graviton Energy 

For the shift in excited energy levels, we need to identify the non-leading exponential behav- 
iors in the free closed string propagator. Since we also want to pick out a specific spin, it is 
iimportant to work with the complete amplitude, including all — 2 factors: 



{N+l,{q^}\0,{q^})tZf = [P'='-'^(iV + l)]^-%'<tr'' 



(30) 



^-T/f/ closed 
^•^•^ total 



To 
' 2 



(9o,Ar+i Qo,o 



\2 M-1 



N + 



+ smh Xl, ( (ql^^^, + ql^,) coth(iV + 1) A 

TO=1 ^ 

29TO,Ar+l ■ Qmfi 



smh{N + 1)A: 



(31) 
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where q denotes a L> — 2 dimensional vector. Here we identify the m < M/2 modes with 
cosine modes and the m > M/2 modes with sine modes. The graviton state on the lattice 
has energy 2X1 above the ground state energy, and involves only the 1 modes. Also since it 
is a symmetric traceless 0{D — 2) tensor, its contribution to the closed string propagator 
resides in the second order term in the expansion of 



exp 



' To{qlN+i ■ Qjo + QIn+1 ' Qjo) smh 
sinh(A^ + 1)A^ 

1 + 2To(g^,^+i ■ + QtA^+i ■ Qio) ^inh AJe-(^+i)^? (32) 
+2T^iQw ■ + Qw ■ QUf A^e-2(^+^)^^ 

as A?" + 1 — > oo. The first term (the 1) propagates an 0{D — 2) scalar, the second term 
propagates an 0{D — 2) vector, and the third term propagates a combination of a traceless 
symmetric tensor, and antisymmetric tensor and a scalar. Since all closed string states must 
be cyclically invariant, the second (vector) term is projected out of the spectrum. But, in 
any case, for the shift in the graviton energy, we may simply drop the first two terms, and 
keep only the symmetric traceless, cyclically invariant part of the third term. 

To identify the contribution of the cyclically symmetric states to the third term we 
consider the new coordinates 



q\ ± iq\ 



(33) 



which acquire the factor e^^'^*/*^ under a cyclic transformation of one step. To make a 
cyclically invariant combination, we must have equal numbers of + and — factors: 



9i,Ar+i ■ 9i,o + 9i,Ar+i ■ 9i,0 ~ 



qVQi.AT+I ■ 9l,0 + 9l,Ar+l ■ 9l,o) 



(34) 



where the last line shows the only contribution that survives the cyclic symmetry require- 
ments. To see the SO{D — 2) content write 



ck „sl 



(35) 



The third (imaginary) term gives the contribution of the anti-symmetric tensor, whereas 
the first two (real) terms give a symmetric tensor, which can further be decomposed into a 
traceless symmetric tensor and a scalar: 



+ 



Jkl 



D-2 



Co + qlo) (36) 
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The quantity in square brackets on the right represents the contribution of the graviton in 
all spin configurations. To identify the graviton energy shift it is sufficient to simply pick 
one polarization with k ^ I so that the trace subtraction drops out. 



exp 



eye mv 



smh{N + l)Af 

1 + T^QtUi^^:N^i ^T,oQti sinh^ A^e-(^+^)^^ + • • • (37) 

= 1 + TM'N+iQtN+i + ?fiv+i??iv+i)(??o9?o + Qto<o) smh^ A^e-^(^+^)^J (38) 

+7^0 (?^5v+i9? Wi - gfiv+i^iVi) (?^o??o - llfilto) sinh^ A?e-^(^+^)^5 + • • • (39) 

The first 1 term can be dropped in the calculation of the graviton and also the anti-symmetric 
tensor energy shifts, since it contributes only for scalar states (the tachyonic ground state 
and the massless dilaton. Then we can take J, L — > oo in 



{N + h{xf}\0,{x^})fZf-\K,J) 

~ Sinh^ A^ ^V^^osed^ j^jyopen^j^^jydosed^^^-^ - ^ 

/ dqldqf[q%qf, + qf^qf^MA + €^<MA,i + <i<i][€< + 

giW'total-(Jf-l)p-2)Bo-Tb[(a!i:,-a=f')2+(a;^-a=f)2]/4 ^^^-^ 

The sum over K of this expression should be compared to the free closed string propagator 
for the graviton 

[V{N + 1)]^-^ sinh^ A^e-2(^+i)^? [q%qf, + q%qf,] [ggg^^ + gggj] (41) 
to read off the graviton energy shift: 



MT^ sinh^ A^ J2 



I 



K=l 



g2KAf-(K-l)(D-2)Bo 



V{N + 1) 

dqldq^\qt.€^ + €i€M1^4^ + gj;g|^;,]e^^---^°[(-M-="i )^+(="M-f 

Ob 

-l^ihk'h' + kl'5lk')^PGraviton + ^hlk't' (42) 



where the C term contributes to the dilaton energy shift. (If C = the graviton and dilaton 
remain degenerate.) Here iW is the same expression (15) that we used in the evaluation of 
the ground state energy shift. A simple way to isolate the graviton shift is to simply choose 
index values for which the C term decouples. For example, take k = k' = 1 and I = I' = 2: 



-(^^Poraviton = 2MT2 siuh^ A^ 



/ 



K=l 



D-2 



V{N + 1) 



dq^dq^[ql\qti + + g^^,,g^^Je'^-'-^°[(-M-f)^+(----^^)^]/^ 

13 



2MT^ sinh^ ^ 



K=l 



q^closed (^J^ q^open ^ q^closed ^j^^ 



D-2 



^2KXl-{K-l){D-2)Bo 



V{N+1) 

J dq^dq^JqUli€^qti + ggi^lc^.i^i^f ije^'^— ^"[(^^-f )^+(^^-f 



2MTj sinh^ ^ 
ft:=i 

,cl „cl \2 , /„sl „sl \2 



D-2 



q^closed (^J^ q^open ^ q^closed ^j^^ 

P(7VTl) 



^2KAf-(A:-l)Bo 



(43) 



where in the second and third forms we take advantage of the fact that the integration is 
over independent Gaussians, so the language of correlations reflected in the (■ ■ ■) notation is 
appropriate. The correlator of cosine modes is just that of the free closed string and is easily 
shown to be 



sinh LXl sinh JAJ 



To sinh A^ sinh(ii: + J + L)A5 2rosinhA5 



(44) 



in the limit J, L ^ oo. 

The correlator of sine modes is of course more complicated because they involve the 
nontrivial overlap of the closed and open string modes. 



RecaU that iW" = -nf 
to evaluate 



A B 



TT 



To sinh Af 

/ dqK,mdqL,mqL,iQK,ie'^" 



(45) 
(46) 



B A 



J dqK,mdqL,me^^" 
q so adding a source term J^q, we complete the square 



A B 
B A 



J 



With the definition 

A B 
B A 



A' 
B' 



we then have 



A' B' 
B' A' 

{A-BA-^B)-^ 
(B-AB-^A)-^ 



-{{A + B)-' + {A-B)-') 
l{{A + B)-'-{A-B)-'), 



2To sinh A^ 



(47) 



(48) 



(49) 



(50) 
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K/M 

Figure 4: Log-linear plot for the summand of the rescaled ground state energy shift 
(51), with each curve exhibiting its dependence on K for a fixed value of M. We 
have rescaled the horizontal axis to K/M in order to demonstrate that for a large 
enough value of the latter, logSP^ develops the same slope. The plot also justifies 
our choice Kmax = 2M for the cutoff in the K sum. 



3 Numerical Analysis 

3.1 Closed String Ground State 

Here we will perform a numerical study of the 1-loop shift in the ground state energy for the 
closed string in D = 26 dimensions, with the help of Mathematica. In particular, we will 
rescale the energy shift (26) by an overall factor of M (notice also APJ^^^^^^ and 5P]^ differ 
by a minus sign). 



aAP, 



G, closed 



M 



E 

K=l 



K=l 



coth(Msinh"^ 1) 
MV2 



1/4 ^yA 



:J{X'f,-X'^)/2-{K-l)Bo 



n 



A/-1 
m=l,odd 



det-1/2 



A 
B 



B 
A 
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(51) 



where A and B are given by (27), (28), as it will turn out that this yields a finite quantity for 
M — )■ oo, which is just what is expected from the bulk term in AP^ . It is instructive to start 
by investigating the dependence of the summand 6P^ on M and K, as a means to also set a 
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Figure 5: Summand for the rescaled ground state energy shift as a function of M 
for K = 2, 3, 4, 5 <C M, including fits of the form (52). 



reasonable cutoff Kmax > K iw the sum. For fixed M, 5P^ decreases rapidly for increasing 
and for ^ M it becomes proportional to e~^-^^®^/*^ times an M-dependent factor. 

This fact is evident in figure 4, which also indicates that Kmax = 2M is a sufficiently large 

cutoff, given that the largest term in the sum is 6Pi = 1 for any M, and SPxmax ^ 10~^^ 

for M > 195, precisely indicating the accuracy of our cutoff. 
For fixed K <^ M, we also find a good fit to 



5P, 



K 



(52) 



where roughly cf ~ K~ and ~ K . The particular examples K = 2,3,4,5 are 
presented in figure 5. Since this behavior changes as K becomes comparable to M, it 
suggests that the sum over K should give rise to a log M/M^ term due to 



M 



^ cf ~ H{M) ~ log M for M > 1 



(53) 



K=l 



where H{M) is the M-th harmonic number. Similarly the sum of cf yields harmonic numbers 
of order 3, whose large M expansion suggests the absence of an 1/M term. 
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We should also note that in both regimes we examined, the general structure of the 
dependence on Af and K is correctly captured by the discretized version of the continuum 
amplitude (9). Denoting the summand of the latter (up to a proportionality factor) with a 
prime in order to avoid confusion, we have for D = 

^""-'^[MsinHnK/M)) = \ ^ _ ^ ^ o(^) K«M ^''^ 

For K ^ M there is also approximate agreement in the value of the exponent, although 
for the K <^ M expansions more detailed comparison of the coefficients of 5P^' and 5P^ 
reveals that they are not simply proportional to each other. 

Armed with this intuition, we proceed to the numerical calculation of (51), summed up to 
Kmax — 2M, and for values of M ranging from 5 to 995 in steps of 10. We fit the generated 
data for different subintervals between M G [195, 995] to ensure that M is sufficiently large 
and to test the stability of our fits, and also calculate the value of as an estimate of their 
goodness. We find that indeed the fit 

c +c ^ +c^°S^ f55l 



M M2 ^ M2 

with 

ci = 1.158863267±3- 10~^ C2 = 2.799 ± 0.011 , cg = -2.800 ± 0.002 (56) 

matches excellently with the data, with the values of the coefficients varying only mildly 
when fitting different subintervals in M (the error estimates are precisely taking this interval 
dependence into account). 

Our main finding of this section, (55), is plotted against the numerical data in figure 
6. We've also included the fit with C3 = to show its insufficiency in accurately describing 
the data. As far as the fit with C2 = is concerned, it leads to values of C3 which may 
differ Tip to 7% depending on the interval of the fit, and generally one should also expect a 
constant multiplying M inside the logarithm. As additional evidence that the A/- dependence 
is indeed correctly captured by (55), we also mention that when fitting the entire interval 
M e [195, 995], the value of B? differs from 1 by a mere 3 • 10~^^, whereas for the C2 = 
and C3 = cases the differences are 6 • 10~^ and 0.002 respectively. Finally, our expectations 
for the absence of an 1/Af term is confirmed by the fact that its inclusion yields unnaturally 
small values for its coefficient, and does not substantially improve the fit. 

As described in the introduction, the bulk counterterm can be chosen to cancel the 
contribution to P~ proportional to and what is left gives the physically significant con- 
tribution. Lorentz invariance requires that this residuum behave at large M as 1/A^, since 
Am^ = 2AfaTo(AP~ — Bulk Term). Our results (55), (56) contradict this requirement be- 
cause of the In M dependence. Taken literally, the result implies a logarithmically divergent 
self mass: InA^ = lnP+/aTo = ln(l/a) -|- ln(P"'"/To), which is to be expected from the dq/q 
behavior in the covariant expression for the self-energy. As is well known this divergence can 
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Figure 6: Rescaled ground state energy shift as a function of M. We present fits 
with and without a log M term, in order to demonstrate the necessity of the latter 
for agreement with the data. 



be absorbed in a renormalization of the Regge slope parameter a' = 1/2'kTq. But the InP"*" 
signifies a noncovariant finite part. 

The origin of this In M factor can be traced to the sum over K of the 1/ K dependence 
we saw in the summand when K « M. The lightcone lattice has cutoff the logarithmic 
UV divergence (small K), but the large K behavior is cut off at i^T = (9(M), because the 
level spacing is of order 1/M, so X^i^ll/-^) = 0{hiM). As mentioned in the introduction, 
the presence of the Bq counterterm offers a way to interpret the lattice calculation that 
avoids the difficulty. By adding a small positive constant Bq ^ Bq + e the cutoff on the 
K sum becomes 1/e instead of M, and the residuum will behave as (l/M)ln(l/e) which is 
still divergent as e — ?■ 0, but remains compatible with Lorentz invariance. Then the ln(l/e) 
can be absorbed in a redefinition of Tq before taking M — )■ oo. As we shall see in the next 
subsection, this same interpretation leads to a zero self-energy for the graviton state. 

3.2 Graviton 

We proceed to investigate how the lightcone lattice handles nontachyonic states by looking at 
the 1-loop mass shift of the spin-2 excitation of the closed string, representing the graviton. 
Similarly to the ground state, with the help of (44), (50), we may rewrite (43) in the J, L — )■ oo 
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Figure 7: Coefficient Cq in tlie rescaled graviton energy shift (57) as a function of 
K/M^ with each curve corresponding to varying K for a fixed value of M. On the 
left hand side we plot the entire range of summation for K/M, it is evident that 
curves for different M values are indistinguishable. On the right hand side we zoom 
into a small region in order to see the mild dependence of the log Cq offset on M. 



"ij 1.0010 - 



Numerical data 
25.132 

M- 





Numerical data 




69.1816 





Figure 8: Coefficient Cq in the rescaled graviton energy shift (57) as a function M 
for fixed = 2, 3 < M. 
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Figure 9: Cq5P^ as a function M for fixed K = 2,3 M . Comparing with figure 
5, we notice that the constant terms are equal, and the coefficients roughly 

opposite. 
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limit as 



aAP, 



Graviton 



M 



^ OO _ CX) 



1 + C^)6P^ 



(57) 



K=l 



K=l 



where 5P^ is the summand of the rescaled ground state shift, defined in (51). 

As 6P^ has been determined in the previous section, the only additional numerical com- 
putation that has to be done is for the coefficient Cq . A preliminary analysis shows that 
for fixed M and varying K this is a rapidly increasing function which for K ^ M becomes 
proportional to roughly e^-^^^/*^, as can be seen in figure 7. However given the behavior of 
6Pj^ in the same regime, their product is guaranteed to converge, albeit more slowly. As far 
as the regime ^ M is concerned, we observe that Cq = 1 for any M, and more generally 
Cq = 1 + c/M^ + 0{1/M^). The first two nontrivial examples K = 2,3 are plotted in figure 
8. The fact that the constant term is independent of K and equal to one guarantees that 
the leading divergence for the ground state and the graviton is the same, as it should for all 
states. 
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Figure 10: ^g^^k ^ ^ function of M. Again fits with and without a logM 
term are presented, so as to demonstrate its necessity for matching with the data. 



Before evaluating the entire sum (57), it is again useful to examine Cq6P^ for fixed K. 
As can be seen in figure 9, for individual K <^ M the latter has an expansion in M of the 
form (52), where cf^ are roughly equal and cf' are roughly opposite between the ground state 
and the graviton. Then the sum in K is depicted in figure 10 similarly described by a fit 
of the form (55), where again the coefficients Ci and C3 are found to be equal and opposite 
respectively within our margins of error, however C2 = —1.93 ± 0.04. This is a first hint that 
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Figure 11: Entire summand of rescaled graviton energy shift (57) as a function M 
for fixed K = 2, 3 < M. 



although additional cancellations occur for the graviton, which may remove the divergent 
logM terms, the lattice regularization, in the absence of the e prescription, still leaves an 
unphysical finite mass shift for the graviton. Moving now to the entire summand (57), from 
figure 11 we infer it behaves as 

i(l + C^05P^ = cf + ^ + O(^) (58) 
with cf ^ cf ~ and cf ~ which is consistent with 

M 



~ M72 for M > 1 , (59) 



E 

K=l 

introducing an additional contribution that changes the C2 coefficient for the graviton, com- 
pared to the ground state. Finally, we find that the fit for the entire rescaled energy shift 
(see Fig. 12) 

M -'^ + '^M^' 

with 

ci = 1.158863276 ± 1.5 ■ 10"® C2 = 0.454 ± 0.004 , (61) 

is in very good agreement with the numerical data^ and our previous quantitative and qual- 
itative observations. 

Our results indicate the absence of the InM Lorentz violating effect we found for the 
tachyon'', but still the K sum with a cutoff of (9(M), leads to the undesirable conclusion 
that the graviton would gain a (necessarily finite) nonzero mass at one loop order. This is a 



^In particular when fitting on the range M G [195, 995], differs from 1 by 10~^. Including an additional 
log M/M'^ term yields an unnaturally small coefficient and does not improve significantly. 

In M divergence in the graviton self mass could not be absorbed in the Regge slope parameter because 
at zeroth order the graviton is massless. 
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Figure 12: Rescaled graviton energy shift as a function of M. 



(more subtle) violation of Lorentz invariance. However, with the e prescription introduced to 
interpret the tachyon mass shift, this difficulty is avoided. Putting Bq ^ Bq + e cuts off the 
K sum at ~ 1/e, so the large M expansion at fixed e encounters no contribution and 

hence no shift in the graviton mass. It is important to appreciate that this interpretation 
requires taking the continuum limit before taking e — ?■ 0. 



4 D-branes 

We extend the discussion to the case when several of the transverse open string coordinates 
satisfy Dirichlet conditions. In current popular terminology this is known as closed string 
theory in the presence of D-branes [13]. To avoid confusion we will call such coordinates 
yj. We shall follow [14] in adapting the lightcone lattice to Dirichlet boundary conditions. 
Starting from the closed string potential energy for one such coordinate 

J, M 

Klosed = -^^(Vi+^-Vif (62) 
i=l 

we pass to the potential energy for a Dirichlet open string with say ?/; = by the following 
substitution 

{yi+i-yif + {yi-yi.if ^ yf^^+yt^ + 2Kyl (63) 
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In other words we keep all M degrees of freedom on the lattice. Instead of trying to set 
= 0, we decouple it from the other coordinates and give it a potential T^nyf that makes 
it produce an energy of 0{1) in lattice units. This means that excitations of yi will have 
infinite energy in the continuum limit, which therefore locks this degree of freedom in its 
ground state. The normal modes for the Dirichlet open string coordinates g^m are defined in 
Appendix A. In this case it is convenient to alter the corresponding normal mode expansion 
for the closed string, which for M odd are 

■{M-l)/2 



1 



Vk 



Qo + 




M 



E 

m=l 



Qcm COS 



2mnk 



M 



+ Qsm sin 



2mnk 



M 



(64) 



The modification for M even can be found in Appendix A. 

We can express the open string normal modes in terms of the closed string ones 

for m even 



Qsm/2 

(M-l)/2 



D 

mm' 



for m odd 



(65) 



m'=0 



where we have defined ^co = Qo/V^, and U^^, is given in Appendix D. 

In constructing the one loop diagram, wc would like the j — 0, N + 1 sites of the open 
string propagator be assigned half the closed string potential energy. We have 



-ToVMiVi + Vm-i + (/« - l)y 



M 



-2U{y). 



(66) 



Thus the loop integrand is given by the product of the three propagators times the factor 
e^(^^ at each vertex. In terms of the closed string normal modes 



Vm 



yi + Vm-i 



M 
2 



qo + 



go + 2 



(M-l)/2 

M ^ 




Q.cm' 

m'=l 
(M-l)/2 



■(M-l)/2 

= V 17 X] 

m'=0 



(67) 



cm' COS ■ 



(M-l)/2 

gem' cos 

m'=0 

(M-l)/2 




m'=l 

2m'7r 



2m'7r 



M 



u{y) 



Qcm' Qcm" 



K 



2m'7r 2m"7i 

1 + cos — — h COS — — — 



(68) 



(69) 



m',m"=0 

where we defined Qco = qo/V2. Then the one loop correction to the closed string propagator 
is 

E / ^yf clyfe-(^-)^o°+^(.-)+^(.^)(L, {/}|0, {y'}r-'{K, {/}|0, {y^})^ 
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(J,{y^}|o,OT) 



K 

where as before we display only one faetor for the Dirichlet coordinate y. Note that we 
expect ^ Bq because of the different boundary conditions. 

4.1 Closed String Tachyon Scattering off D-brane 

For the case where yi — yf — we have, for each Dirichlet coordinate y, a term 



iW + U(y'') + U(y'-) = 



+ J^^Qc^ + ^«-) ^^^^ -^^^ + J2^1cm + qZ) sinh coth LA^ 

+ ^{QDm + 9Dm) Sinh A° coth KX^ - 2 ^Om^Dm- ^"^^ 



sinhi^AD 

m m '"■ 

2 ^^^^"^ t K K T T \ f 2m'7r 2m" n\ 

2^ [Qcm'Qcm" + 5cm'?cm" ) U - 1 + COS + COS 1 

m',m"=0 



The next step is to change integration variables to the closed string normal modes, Qq — 
qco\/2, qcm, Ism- The Jacobian for the change of variables yk qo, qcm, qsm is unity, and 
further changing go ?co gives a factor y/2. The equality qD2m — qsm means that integrating 
over the closed string sine modes simply implements closure on these modes. Thus we can 
write 

{N + l,{xf}\0,{x^})tZ' 

= pciosed(^ ^ 1)^^^^^^^)^^^ 1 2rf,lci,le^^'''-(--)-o° (70) 



where 

iW 



D/ ^ 

2 



Qo^j + Qo^i + 'icm sinh A^ coth JA^ + ^ g,^^ sinh A^ coth LAJ 

m m 

M M 



+ E (^DL + ^D^)«-hAj^cothi^Aj^-2 J] q^^q^ 



sinhXAD 

m=l,oad m=l,odd 

2 ^■'^^J^/^ t K K L L \ f 2m! Ti 2m" 'k\ 1 

~ M ^ {^cm'Qcm" + Qcm'Qcm") U - 1 + COS -j^ + COS -j^ 1 (71) 
m',m"=0 ^ ^ - 
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Taking J, L large, the factors in front of the integral reduce to 



N 



\/2, 



n 



m,oddSinh AD 



nLlT^^/^sinhA: 



(M-l)/2 

^ n -«-hA: 

27r TT 

m=l 



(72) 



Meanwhile 



(M-l)/2 

J] fem' + 9cm)smhA; 

m'=l 

(Af-l)/2 



To 
" 2 



m' m"=0 



— — K — 1 + COS 



M 



+ 5^ ilcm'Qcm" + QL'Qcm") ^mm'^^mm" Smh A° COth /TAJ^ 



m=l,odd 



M 



2m'7r 



+ cos 



2m"7r 



(M-l)/2 



M 



M2 



V «^ «^ V r/D /"/D 

/ ^ Hcm'^cm" / j ^ mm'^ mm" 



sinh A^ 
sinhXAE 



(73) 



m',m"=0 m=l,odd 

Just as in the Neumann case it is convenient to absorb the factors in square brackets into a 



rescaling of the integration variables qm — <lcm\/ (^o/tt) sinh A^ for m = 1, • • • (M — l)/2 and 

/2(Ar + l)e^(^-'i^'^^' ^5n-E^.odd A° )/2-(K-i)B° 



go = qQ\/TQ/2Ti = qcoy/To/TT. 
{N + l,{xf}\0,{x'})f,Z^' 



pclosed(7Y_^ 1) 



JL 



rim.odd sinh AD 



AD 



(74) 



- sinh ^ 1 — sinh ^ ' ^ 
2 



In 



1 + ^2 



2 2 1 + « + y«(«; + 2) 



(75) 



Comparing the prefactors in this formula with the corresponding factors for the Neumann 
case we see that the extra factors are 



2{N + 1) e-^^M/2^sinh A^ 



JL 



N+1 
JL 



' 1 - 772^ 



(76) 



where for brevity we have defined rj = 1 + k — a/k(2 + k) ^ 2 - 73 ^ 0.268 for k = 1. 
The •\/(A'" + 1)/JL factor just reflects the fact that the intermediate open string has its ends 
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fixed in space. In the Neumann case this factor would instead be 1. The factor rj^~^ can be 
absorbed in the boundary countcrtcrm, converting its zero couphng value back to Bq. The 
expansion in powers of rj represents excitations of order 0{1) in lattice units, which will be 
suppressed in the continuum physics. 

Finally we turn to the matrix determinant that results from the execution of the Gaussian 
integration. For this we need to spell out W^' which remains after integrating out the closed 
string sine modes. Expressed in terms of the new variables q^/ , we write 



— TT 



rn,m" 

M 



^00 



4° — 



M2 



m=l,odd 
M 



M 



2 ^ j3 j3 sinhA^coth/TA^ K + cos2m'7r/M 
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M 



""^ ^sinh A^, My^sinh A^, 

sinhAHcothXAS 



(77) 
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2 ^^i:^dd -"V^nsA^V^ShA^ 

K - 1 + cos(2mV/M) + cos(2m'V/M) 
My^sinh A^,y'sinhA^„ 

2 V" r^D rrD sinh A° 



-°0m' — -°m'0 



M2 
2 

2 



m=l,odd 
M 



'sinhXAS 



V r7° r/^ - 



sinh A!S! 
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V - 

/ J mm' mm" . 



sinh i^Aj^-^sinh A^, 
sinh AE 



sinh XA^^sinh A^, ^sinh A^„ 



(80) 
(81) 

(82) 

(83) 



m=l,odd 

For a Dp-brane there are D — p — 1 — > 25— p coordinates satisfying Dirichlet boundary 
conditions. Putting everything together we have for the zero energy amplitude for a closed 
string tachyon scattering off a Dp-brane: 



— aMa, closed 
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1 + K — k{k — 1) 



cothMsinh-^ 1\ e-^^^^=i 



Ma/2 

A B 



n 



m=l,odd 



24 



- p-1 






\l MTo'\l 



^det-V^M" 
2^ aet I 



25 -p 

(84) 
(85) 



where the scattering amplitude is obtained from the one loop correction to the two closed 
string function by stripping off the factor ^y MTq {N + 1)/ 2nJL for each Dirichlet dimension, 
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as explained at the end of Appendix C. The K = 1 term is not included in the scattering 
amplitude since it contributes to the / term of the S-matrix. 



4.2 Graviton Scattering off D-brane 

Let us take the graviton polarizations to lie within the Dp-brane. Then in parallel to the 
derivation of (43) we must simply insert the factors 

2T^ sinh^ mqt,qt^f + {<lUtr?)e'^''^ = ^(1 + {B[,e^'^') (86) 
into the K summand for the closed string tachyon scattering amplitude. 



5 Numerics of D-branes 

In the case of a closed string tachyon scattering off a Dp-brane, it is convenient to define the 
quantity 

/ dei{A + B) dct{A - B) 
y Mdet(A^ + fi^)det(A^-B^) ^ ^ ' 

such that the corresponding zero energy amplitude (85) may be rewritten as 



aM 



G,closed 



M 



K=2 




25-p 



T.^PK(\-r—k'-) ' (88) 



where 6P^ is the summand of the tachyon energy shift (51). For simplicity we have set 
To = 2n and also k = 1, as we have checked that varying its value does not substantially 
change our results. 

Hence the ratio essentially encodes the difference between Neumann and Dirichlet 
boundary conditions, and we find that for ii' 3> M it falls off exponentially to a value of 
(9(1) which depends very mildly on M, as can be seen in figure 13. For fixed K <^ M we 
find that the fit 

matches very well with the data (see figure 14 for indicative values of K), and by further 
examining the fitted coefficients for different values of K, we infer that they roughly vary as 

K 0.115 , , 

of ^ 0.724 + — — -, (90) 
K — 1 

-0.204 + 0.369i^, (91) 



see also figure 15. By factoring out the constant term in of-, which dominates rx, we 
can infer that the leading dependence of the entire amplitude (88) on p will be equal to 
(0.724^1 - 7^2)25-? or equivalently 0m7^^-P. 
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Figure 13: Linear plot and Log-linear plot of tk and its first difference respectively, 
where each curve corresponds to fixed M and varying K. On the left we see that 
for K ^ M rx approaches a value which only depends on M, and on the right in 
particular that the difference of rx from this value is proportional to g-6.28ft:/Af 
any M in this regime. 
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Figure 14: as a function of M ioi K = 2, 15, including fits of the form (89). 
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Figure 15: Values of the fitted coefficients af-, af- in (89) for K = 1, 
against their estimates (90)-(91). 
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Figure 16: Rescaled zero energy amplitude for a closed string tachyon scattering 
off a Dp-brane for p = 3, 10, 15, 20 as a function of M. 



In fact, it is evident that the only 1/{KM'^) term in the summand (88) can arise when 
the aforementioned p-dependent factor multiplies the (9(1/M^) term in 6P~ (52), and given 
(53), our analysis constitutes a definite prediction for the relation of the log M/M^ term of 
the entire sum in the absence or presence of D-branes. Furthermore, the structure of is 
such that no l/K"^ or K'^^~'^ /M"^^ terms appear in the summand (88)®, and consequently no 
1/M terms appear in the sum, so that its expansion to next-to-next-to- leading order in M is 
expected to be of the the same form as in (55). Having gained this insight from the analysis 
of Tx and the summand, we proceed to evaluate and fit the entire zero energy amplitude 
(88), and find that its leading large M behavior is indeed captured by an expansion of the 
form 

aMcdosed , 1 logM 

I^='' + ''M^^''^- ^'^^ 

In figure 16 we compare the fit with the numerics and for sample values of p, and give 
for these cases the values of the p-dependent coefficients. It's worth taking a closer look 
at the dependence of C3 on p in order to see the effect of the presence of D-branes on the 
undesirable logarithmic divergence, and also verify the prediction for its value based on the 
aforementioned rx analysis. Plotting log(— C3) as a function of p G [1, 25] (see figure 17) we 



would require that af^ has at least a linear term in K, whereas the K^/A/P"' terms that appear 
always have I < n — 3 for any n. 
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Figure 17: Log-linear plot of the log M/M^ coefficient C3 in the large M expansion 
of the tachyon-Dp-brane scattering amplitude (92), as a function of p. 



identify a clear linear dependence, which we can fit in order to obtain 



This is precisely the coefficient of the corresponding term in the closed string tachyon energy 
shift (56) times the factor we identified below (91). The dependence of C3 on p shows that the 
presence of D-branes only softens the divergence in the sense of reducing its coefficient, but 
cannot remove it completely. Furthermore, these results serve as additional evidence that 
the existence of the divergence is solely due to the summation of the term in the summand 
of the closed string tachyon self-energy, which behaves as 1/i^' for small K. 

Finally, let us conclude by briefly examining what changes when instead of a tachyon 
we have a graviton whose polarizations lie within the D-brane. We flrst recall that in the 
absence of D-branes, the summand for the graviton rescaled energy shift (58) has no 1/M^ 
term, meaning that the graviton remains massless at short distances. On the contrary here 
we notice that the presence of such a term in tk (89) will carry through to the graviton 
summand, which can in turn be interpreted as mass generation due to the explicit breaking 
of Lorentz invariance by the Dirichlet boundary conditions. Although in this particular case 
the mass is tachyonic, generally having a mechanism for mass generation may be viewed a 
desirable feature, as we are ultimately interested in using open string theory to probe QCD 
phenomena. Thus the massless spin-1 states of the former will have to acquire a mass if they 
are to be put in correspondence with massive gluonic states. 



C3 = -2.800 ■ (0.697) 



25-p 



(93) 
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Figure 18: Ratio of the first two coefficients in the large M expansion of the 
summand of the tachyon-Dp-brane scattering amphtude (95), as a function of 
p. Each hne represents the i^'-th term in the sum, K = 2, . . . , 10, and the 
value of the corresponding slope is also provided in the legend. 



In more detail, if our previous fits (58)- (89) have captured the M-dependence of each 
factor correctly, we should have 



aM 



Graviton 



M ^'2 

K=2 



E[^(l + ^c)5^^]( J^^-x) (94) 



OO 



K=2 

oo 



25-p 



= E + (25 - P)cf (afT'-'of -rij + 0(-^)] • (95) 

K=2 

In particular the ratio of the second to first term should be a linear function of p with slope 
of /af^ which always becomes zero for p = 25, as we will now explicitly confirm. We fit the 
summand in (94) for p = 1 . . . , 25 and i^' = 2, . . . , 10 to inverse powers of M up to (9(1/M^), 
take the ratio of the coefficients of the first two terms and perform linear regression in terms 
of p for each value of K. As can be seen in figure 18, the dependence of the ratio of the 
coefficients on p is clearly linear with the right intercept, and comparing the slope to a'^ /af 
we find that in all cases it differs less than 0.1%. 
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6 Discussion and Conclusion 



In this article we have made a modest beginning to a critical study of the effectiveness of 
the worldsheet lattice, introduced in [3], for implementing a regulated bosonic string loop 
expansion. We have limited our analysis to the one-loop corrections to the closed string two 
point function in the presence of D-branes, but we plan to analyze the open string self-energy 
in a second paper. 

A convenient way to summarize our results is to recapitulate the lightcone lattice in- 
terpretation of the divergences in the well-known covariant expression for the closed string 
tachyon self-energy (6): 



The first term dq/q^ gives a quadratic divergence, which in the covariant description is 
associated with the closed string tachyon disappearing into the vacuum. On the lightcone 
worldsheet lattice, we have seen that this divergence is interpreted as a contribution to 
P~ ~ — ckM, which can be cancelled by the bulk counterterm proportional to the area 
M{N -\- 1) of the worldsheet lattice. The only states that survive the continuum hmit are 
those with the smallest value of a. Since only energy differences are physically significant, 
this shows that the quadratic divergence is physically inconsequential. 

The second term dq/q gives a logarithmic divergence, associated with the disappear- 
ance of a closed string dilaton into the vacuum. This divergence can be absorbed into 
a renormalization of the Regge slope parameter a' . Our analysis using the worldsheet 
lattice has shown that this divergence shows up as a contribution to P~ behaving as 
(aM)^^ InM = (Tq/P^) \n{P'^ / aTo) . The noncovariance of the finite residuum was caused 
by a cutoff on the K = T/a sum proportional to M. If the cutoff were M independent the 
noncovariance would disappear. We then noted that such a cutoff is naturally introduced 
by adding a small positive number to the boundary counterterm Bq ^ Bq + e. Then the 
InM is replaced by ln(l/e) which can then be covariantly absorbed in a renormalization of 
a', before the continuum limit M — )■ oo. 

This same e prescription also prevents the graviton from gaining a mass in the absence 
of D-branes. Interestingly, in the presence of a D-brane the low energy amplitude for a 
graviton scattering off the D-brane, which can be thought of as a "self-energy" for a graviton 
propagating parallel to the D-branc. suggests a nonzero graviton "mass" in spite of the e 
prescription. This is a consistent outcome since Lorentz invariance is broken by the Dirichlet 
boundary conditions. This is important in applying these ideas to large N QCD, since the 
closed string is supposed to model glueballs, all of which should be massive. 

At a more fundamental level, we can recognize e as a natural and bona fide physical 
parameter of the lightcone worldsheet system: it is simply a measure of the boundary energy 
B = Bq + e associated with the disappearance of a link. It is also the minimum energy 
assigned to a free string bit. The free open string has a Lorentz invariant spectrum for only 
one value B = Bq (or e = 0), but it makes sense to study the physics of the system as a 




(96) 
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function of B (or e). For sufficiently large B, 





(97) 



M free string bits will have an energy greater than the ground state energy of a closed string 
of size M. In fact, for B this big, any system of open strings with total bit number M has 
energy larger than -PcTosed.G- Thus for S > 2G/7r, the free closed string is stable against decay 
into any number of open strings, and our scheme to sum planar diagrams should make good 
physical sense. 

Dialing the value of B gives us a new tool to analyze the fate of open string tachyons 
in the bosonic string on the light cone worldsheet lattice. One can attempt to simulate the 
sum over all patterns of missing links, as explained in the introduction, in the closed string 
propagator for a range of B values, and then study how the physics changes as S is gradually 
decreased. If the dynamics is favorable, the system should find the true vacuum with all 
traces of tachyons removed. Of course it is possible that the dynamics does not stabilize the 
theory lending weight to the prevailing opinion that the presence of tachyons of the bosonic 
string theory is an incurable disease which can only be cured by replacing the bosonic string 
with the supcrstring. 

By introducing a D3-brane and suitable orbifold projections, it is possible to arrange the 
gauge boson sector of the open string to enjoy the same dynamics as gauge field theory in 
4 space-time dimensions. Of course the open and closed string tachyons are still present, 
but they could simply be a symptom that perturbation theory is being attempted about the 
"wrong" vacuum. If so, analyzing the worldsheet lattice system as a function of B could 
provide a way to find the "right" vacuum. It will be an interesting exercise to apply these 
techniques to the problem of quark confinement. 

Acknowledgments: This research was supported in part by the Department of Energy under 
Grant No. DE-FG02-97ER-41029. 

A Normal Modes 

A string with P"*" = MgTq is described at a fixed time by M coordinates Xi oi yi,i = 1, • • • M. 
In this article we require several normal mode decompositions depending on the boundary 
conditions. 

Neumann Open String 



1 



?o + \/ 77 




m=l 



m7r(i - 1/2) 
M 



(98) 




(99) 
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Closed String (Neumann) 
M odd : 



1 



M even : 



m'"'' Vm 




(M-l)/2 

E 

m=l 



2m7r(i-l/2) 2m7r(i-l/2) 

gem COS — h Qsm Sm - 



M 



M 



(100) 



{qo + g5M/2(-)*) 



+ 




M/2-1 

E 

m=l 



2m7r(i-l/2) . 2m7r(i-l/2) 

gem COS ^— — + Qsm Sm - 



M 



M 



(101) 



2 2m7r(^-l/2) A2"v^ 2m7r(^ - 1/2) , ^ , 

— ^ cos , Qsm -\ Tjl^^i Sin ^ (102) 




M 



M 



qsM/2 = V7^X^(~)*^'' 



?0 



i=l 



i=l 



(103) 



Dirichlet Open String 



M-l 



/ 2 . mnk . , ^ ,^ -, 

Vk = \Jj^Z^ ^Dm sm for A; = 1, • • • , M - 1, 

m=l 
■M-l 



M 

2 



QDt. 




M 

mnk 



Vm — Qbm 



M 



Vk sin < m < M, goM = Vm 



k=l 



(104) 
(105) 



Closed String (Dirichlet) 
M odd : 



Vi = 
M even : 

Vi = 



1 n-{M-i)/2 
\ — ^ 



M 



m=l 



gem COS 



2m7ri 



^(go + gcM/2(-)0 + V^ E 



M 
2 



M 

M/2-1 



+ ?sm sin 



2m7ri 



m=l 




M 



cos 



2mm 
1^' 




gem COS 



2mm 



M 



+ g^m sin 



2mTii 



M 

M 



sm 



2mm 



i=l i=l 



(106) 

(107) 
(108) 
(109) 
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B Determinants 

Define the normal mode frequencies of a one dimensional harmonic chain 



4 sin 



9„ = 4 sin^ 



2(Ar + l) 
mTT 



7jfc = 4 sin' 



2M' 
2(^ + l/2)7r 



m = 0, 1,...,M- 1 

A; = 0,l,...,ir-1 



mTT 



4sin^ , m = 0, 1,...,M- 1 

M 



(110) 
(111) 

(112) 
(113) 



where a, ^, 7, 5 are the modes of a Dirichlet-Dirichlet, Neumann-Neumann, Dirichlet- Neumann, 
closed chain respectively. Then we are interested in the following determinants: (see, for ex- 
ample [15]: 



N M-1 



iV M-1 



D 



DDDD 



1/2 



D 



DDDN — 



D 



DD ring 



n nK+^' 

n=l m=l 
M-1 K-\ 

n n(/^-+7^)"'' 

m=l fc=0 
N M-1 



DNDN 



n nK+^-)"'^' 

n=l m=0 
M-1 ftT-l 

m=0 fe=0 
N-\ M-1 



(114) 



ND ring 



n=l m=0 



n=0 m=0 



where the subscripts denote Dirichlet (D) or Neumann (N) boundary conditions on each of 
the four edges of the rectangle. The cylinder determinant with Neumann boundary condi- 
tions must be defined to exclude the overall zero mode: 



M-1 

i^NNring = ^MTV J] (/9„ + U"'^' = H ^'^DD 

(n,m)7^(0,0) m=l 

— Y j^^-^DI^ "ng 

The following product identities can be easily derived: 



ring 
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n sin(A^ + l)«; -p-r 
sm K 



n=l 
M-1 



m=l 



sin2(MK/2) 
sin2(fi;/2) ' 



k=0 
M-1 



m=l 



cos[(2ir + 1)k/2] 
cos[k/2] 

sin Mk. 
sinK 



(116) 
(117) 



where z and k are related by 2; = 4 sin^[K/2]. Applying these identities at z — 0, k — shows 
immediately that -Ddndn = ^dddd/ VN +1 and -Ddnnn = -Ddddn- 
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For convenience we collect here the expressions for the open and closed string propagators 
at vanishing initial and final coordinates [3]. These quantities are the determinants just 
discussed in their various guises. For example 'P°P"=n is simply related to ^^dndn, with one of 
the products performed using the identities (117). 



^oTd"(^ + i) 



1 /T 

* ^ ' m=l 



sinh(Ar + l)A^ 
sinh X° 



-1/2 



rr. \ (Af-l)/4 M-1 

£o_ 

2tt 



n 



m=l.odd 



sinh(iV + 1)A^ 



sinh X" 



1 



Tn 



M/2 M-1 r- . 



n 



m=l 



sinh(A^+ 1)A^ 



-1/2 



.-1/2 



sinh As 



-) 

277 J 



(M-l)/4 (M-l)/2 



n 

m=l 



smh{N + 1)A 



sinh A£, 



-1/2 



^^Z^T^Tm (118) 
(119) 

(120) 

(121) 



'T-Nclosed'T-Nclosed 

— cos sin 



where, for simplicity, we have written these formulas assuming M is odd. If M were even, 
appropriate adjustments to the limits of the products must be made. 



C Propagators 

C.l Neumann Open String Propagator 



{N + l,{xf}\0,{x'}yP^'' = popen^TV + l)e^^°^«" 

^ ' ' open 2 
M-1 



(122) 



N + 1 



+ E + </) -°th(iv + i)a;; - '^-^^f^^) J(i23) 

m=l ^ \ ! m / 



2 sinh ^ ( sin 



m7r\ 



2M 



(124) 



Where the g^'s are the normal mode coordinates for the x's. The right side is the result of 
doing the integrations over all the xl with i = 1, • • • , M and j = 1, - ■ ■ N. The propagator 
spans N+1 time steps and this result corresponds to assigning half the potential energy 

" .3 

i+1 



To - Afl^ to time j = and half to j = + 1. 



C.2 Dirichlet Open String Propagator 

The Dirichlet open string propagator over a time oiK = N + \ steps is evaluated to be 

(Ar + l,{g/}|0,{g*})° = P°(Ar + l)e''^° (125) 
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where 



To 
' 2 



M 



m=l 



Yl («Dm + ?Dm) s^h A° coth KX^ - 2g^^g^, 



sinh A 



m \ ' 



sinhXADy 



126) 



P°(A^ + 1) = 

A^ = 2sinh-^ 



m=l 

2' 



sinh A£ 



AS = A^ =2sinh-^ 



sm 



2M' 



(127) 



■M- 1(128) 



We recall that the above expressions give the the result of integrating over all the variables 
ul, for j = 1,---,N, with half the potential energy assigned to j = 0, + 1, which is 
consistent with the closure requirement. 



C.3 Closed String Propagator 



{N + l,{xf}\0,{x'}) 



closed 



iWr, 



closed 
M-l 



{qoj - qo,iy 



(129) 



To 
' 2 



N+1 



+ E [{qi, + qij) coth(iV + 1)A^ - |l30) 



2 sinh ^ ^sin 



rmr 



(131) 



Where the the normal mode coordinates for the x's. When we divide the closed 

string normal modes into sine and cosine modes, we arbitrarily call the m > M/2 modes 
sine modes and the m < M/2 modes cosine modes. When M is even, the M/2 mode is not 
doubled. The right side is the result of doing the integrations over all the xj with i — 1, - ■ ■ ,M 
and j — 1, - ■ ■ N. The propagator spans iV + 1 time steps and this result corresponds to 
assigning half the potential energy Tq to time j = and half to j = + 1. 

In sums like these it is understood that x-'j^^^^ = x\. Whenever we concatenate at a time j 
propagators with different numbers of missing links, we will understand that we add terms 
To(Ax)^/4 in the exponent so that the potential assigned to time j is that of the system 
with the least number of missing links. For example, the concatenation of an open string 
propagator with a closed string propagator entails the addition of To(x'^ — x\ f/A to the 
exponent. 

Finally, we resolve the zero mode dependence of the propagators in momentum space 



f dp p^T/2P+ i{x'f^-xf^)p ^ fP^ P+(^C 

7 27r V 27iT 



CM ^CM\2 



MTr 



27r(A^ + l) 



^-Mro(a;C'M_j.CM)2/2(Ar+l) 
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f ^ p^T/2P+^{xf'^^-xOM)p ^ /_MIo_ To(5o./-90,,)V2(iV+l) (-^^2] 

J 27r Y27r(A^ + l) ' ^ ^ 

where x^^ = 'Ylk^k/M = qQ/\fM is the center of mass coordinate. From this we see that 
in extracting eigenstate amphtudes from propagators defined with Dirichlet conditions on 
initial and final states, we must not include the factors -^f MTq/2ti{N + 1) . 



D Overlap Formulas 



Neumann Open-Closed 



Umm' 



Qcm/2 



2_ v^(M-l)/2 



m even 



Qsm' Umm' rn odd 



sin(m'7r/M) cos(m7r/2M) 
sin^(m'7r/M) - sin^(m7r/2M) 



(133) 
(134) 



Dirichlet Open-Closed 

1 ry(^-l)/2 ry(M-l)/2 

^^'^ " 7^^° + Vm ^ ^^"^' = Vm ^ ^^^^^ 

^ m'=l m'=0 

gom = qsm/2, for m even (136) 

y/2 '^^-^^ rnvr/c 2 ^ -^--l^ ^-^^ mvrA; 2m'iTk 
QDm — -ttQo / sin — — — h — > Ocm' / sin —-— cos — -- — 
Hum ^ M M ^ ^ ^ M M 

k=l m'=l k=l 

2 ^ '^-l^ . mTr/c 2m' nk /-.^^^ 

= TT > 9cm' / Sm— — — COS — — — (137) 

m'=0 k=l 

where, for convenience, we have defined qco = go/'s/S- The sum over k is easily done 

. mnk 2m'7Tk 5m odd sin(m7r/M) ^ 

2^ sm — — cos — — = — — — — ■ 2/ ; = oddUmm' 138 

MM 2 sm (m7r/2M) — sm (m'TT/M) 

for odd m < M. Wc can unify the treatment of the m — M mode by defining U^^, — \jMj2 
for m' = 0,l,---,(M-l)/2: 

for m even 

^Dm = <( 2 ^"U^/' . (139) 

m'=0 
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XI Icm'^mm' fo^" O^d 



Open-2 Open 



fW I — 2 — 2 

m'=l m'=l 



■M-l ^ M-l 



m'=l m'=l 

Ml 



(1) ^ m'vr /. 1\ mn f . l\ 
I^LL, = > cos \ t cos — \ t 

^ 2/ MiV 2y 

_ sin(m'7rMi/M) sin(m'7r/2M) cos(m7r/2Mi) 

^ sm^{m'n/2M) - sin2(m7r/2Mi) 

i=l+Mi ^ ^ ^ ^ ^ 

_ 1 sin(m'7rMi/M) sin(m'7r/2M) cos(m7r/2M2) 
~ 2 sm2(m'7r/2M) -sin^(m7r/2M2) 



(142) 



(143) 



and we note the identity q^^ \J M\ + ^q^"* \J = qoVM, as expected from the fact that qo/VM 
is the center of momentum of the open string. 

We can also express the g's in terms of the q^^\ gr^^^'s: 



(1) /Ml (2) M2 

^° ^ ^° V M+^° Vm ^^^^^ 



?m' — 




«f/i;4+V2^,g)f/(;<„) (145) 



m=l 



E Robustness of numerical results 

In this appendix, we perform additional tests which verify the correctness of our numerical 
results, within the stated accuracy. In particular, we recalculate dP^ with more significant 
digits^ for indicative values of M and K, and examine the improvement in its accuracy. 

We will use the difference in 6Pj^ between the calculations with different number of 
significant digits, A(5P^), as a measure of our numeric error. We first notice that A{5P^) 
decreases almost exponentially when K increases for fixed M, as can be seen in the example 
of figure 19. 

The decrease is so rapid that we can clearly consider the error of the K = 1 term as 
the error of the entire sum, giving the rescaled ground state energy shift A Pq closed (^l)- 



^In particular, we now keep 20 significant digits compared to the 16 we had initially. 
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M=205 
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1x10-14 
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Figure 19: Log-linear plot of the change in accuracy A{6Pj^) as a function of K for 
M = 205. 
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Figure 20: Plot of the leading change in accuracy A{6Pi) as a function of M. 
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Fortunately, in this case we know the exact value SP^ — 1 for any M, allowing us to also 
obtain the exact deviation of our numerical results from it, which we present in figure 20. 
We first notice that the deviations are centered around zero, implying that no systematic 
error that offsets the value of 6P^ is present. Furthermore, it is evident that the deviation 
increases with M, and for the range M e [195, 995] we have based our fits on, it hes between 
10-11 - 10-10. 

Let us now compare this with the error in the fit for ^Pq dosed (^^) ^ result of the 
uncertainty in the coefficients (56), 

Acl~10-^ a(^)^10-«, a(^^^)~10-« forM~10^ (146) 

where we estimated the smallest possible contribution of the last two terms by replacing M 
with roughly the largest value we used in our numerical analysis. Clearly these uncertainties 
are at least one order of magnitude larger than the errors due to our choice for the number 
of significant digits, so the effect of the latter on the determined values for the coefficients 
Cj will be negligible. This successfully completes the investigation of the robustness of our 
numerical results. 
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